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Abstract 

We study the proposal that the cosmic ray primaries above the Greisen-Zatsepin- 
Kuzmin (GZK) cutoff are gluino-containing hadrons (g-hadrons). We describe the 
interaction of g-hadrons with nucleons in the framework of the Gribov-Regge ap- 
proach using a modified version of the hadronic interaction model QGSJET for 
the generations of Extensive Air Showers (EAS). There are two mass windows 
marginally allowed for gluinos: nrig < 3 GeV and 25 < rrig < 35 GeV. Gluino- 
containing hadrons corresponding to the second window produce EAS very different 
from the observed ones. Light ^-hadrons corresponding to the first gluino window 
produce EAS similar to those initiated by protons, and only future detectors can 
marginally distinguish them. We propose a beam-dump accelerator experiment to 
search for ^-hadrons in this mass window. We emphasize the importance of this 
experiment: it can discover (or exclude) the light gluino and its role as a cosmic ray 
primary at ultra high energies. 

PACS numbers: 98.70. Sa, 14.80.-j 
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1 Introduction 



Since long time light gluinos have attracted attention as a possible carrier of the very high 
energy signal in the universe. In the 80s, they were studied as a possible primary particle 
from Cyg X-3 |l|, 0, now as a primary particle of the observed Ultra High Energy Cosmic 
Rays (UHECR) @, |. 

The observations of UHECR with energies above 10 20 eV impose a serious problem 
(see [f| for recent reviews). 

The data show the presence of a new, nearly isotropic component in the UHECR flux 
above the energy E ~ 10 19 eV ||. Since the arrival directions of the UHECR show no 
correlation with the galactic plane and the galactic magnetic field cannot isotropize par- 
ticles of such energies, this component is thought to be extragalactic. On the other hand, 
the signature of extragalactic protons, the Greisen-Zatsepin-Kuzmin (GZK) cutoff || at 
E ~ 6 ■ 10 19 eV, is not found. The other natural UHE primaries, nuclei and photons, 
must also suffer a similar cutoff. Meanwhile, four different UHECR experiments || do 
not show the presence of such a cutoff. The two highest energy events are detected by 
AGASA |7| and Fly's Eye § at energy 2 ■ 10 20 eV and 3 ■ 10 20 eV, respectively. The total 
number of events with energy higher than 1 • 10 20 eV is about 30, 17 of which are detected 
by AGASA |5J. The accuracy of the energy determination is estimated to be better than 
20-30%. The energies of the two highest energy events |7j and || are determined very 
reliably. 

To resolve this puzzle, it seems that new ideas in astrophysics or particle physics are 
required. 

The proposals involving particle physics include UHE particles from superheavy dark 
matter |ff0l and topological defects p] , the resonant interaction of UHE neutrinos with 



dark matter neutrinos [12], strongly interacting neutrinos [13], new particles as UHE 
primaries [H |4], Q and such a radical possibility as Lorentz invariance violation |Tj| . (For 
more references see also the reviews cited in ||.) 

In this paper, we shall consider a gluino- containing hadron (g-hadron) as a carrier 
of the cosmic UHE signal, being inspired by the correlation between AGN and arrival 
directions of UHE particles suggested by the analyses in Refs. j|, [LlJ. This correlation 
implies that the signal carrier is neutral and is not absorbed by the CMBR. The light 
gluino is a suitable candidate for such a primary: it can be efficiently produced in pp- 
interactions in astrophysical sources, it is not strongly absorbed by CMBR (see below) 
and it produces EAS in the atmosphere very similar to those observed. 

Heavy gluinos are naturally produced in decays of superheavy particles [17| . 



We shall study here the interaction of both light and heavy gluinos with nucleons at 
UHE. In most interesting applications gluinos must be light (see below) . To be a suitable 
primary of UHECR, the g-hadron should satisfy three conditions: 



1. The longitudinal shower profile of Fly's Eye highest energy event with E = 3-10 20 eV 
is well fitted by a proton |L8], |I9|| , though Fly's Eye collaboration does not exclude a 
photon as a primary ||. Therefore, g-hadrons should essentially mimic proton (or 
photon) induced air showers. 
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2. To shift the GZK cutoff to higher energies, the new hadron should have a mass 
in excess of the proton mass: the threshold energy for any energy-loss reaction on 
microwave photons increases with increasing primary mass, while the fraction of 
energy lost per scattering decreases. Moreover, it is desirable that its cross-section 
for interactions with CMBR photons is smaller than the proton's one. This can 
be achieved if, e.g., the mass of the first resonance X that can be excited in the 
reaction g-hadron+7cMBR — ► X is relatively large. 

3. The primary has to be stable or quasi-stable with lifetime r > 10 6 s(m/GeV) (L/Gpc) 
in order to survive its travel from a source (e.g. AGN) at distance L ~ 100-1000 Mpc 
to the Earth. 

In principle gluino-containing hadrons (g-hadrons) could satisfy the above require- 
ments. Below we shall shortly review the status of g-hadrons as UHECR signal carrier. 

To satisfy the third condition, the gluino should be the Lightest Supersymmetric 
Particle (LSP), or have a very small mass difference with the LSP. It also can be the 
second lightest supersymmetric particle, if the LSP is the gravitino; in this case the 
gluino decays gravitationally and its lifetime can be long enough. Theoretically the best 
motivated candidates for the LSP are the neutralino and gravitino. While in minimal 
supergravity models the LSP is the lightest neutralino (in some part of the parameter 
space it is the sneutrino), in models with gauge-mediated SUSY the LSP is normally the 
gravitino. In Farrar's model [ 20| , the gluino is the LSP because the dimension-three SUSY 



breaking terms are set to zero. A theoretically more appealing scenario containing a light 
gluino was developed in Refs. |22|. There, the gluino with mass 1-100 GeV was found 
in a SO(10) model with gauge-mediated SUSY breaking and Higgs-messenger mixing. In 
this model either the gluino or the gravitino is the LSP. In the latter case, the gluino can 
decay but has a sufficiently long lifetime to be a viable UHECR primary, r ~ 100 yr. 

In a physical state, the gluino is bound into colourless hadrons. What is the lightest 
state of gluino-containing hadrons? 

In the 80s (see @), it was argued on the basis of QCD sum rules that the glueballino 
gg is the lightest g-hadron. The lightest baryonic state, gluebarino, was demonstrated 
to be the guud-ha.dion p3| . Gluebarino is a long-lived particle because its decay needs 



either violation of baryon number or R-parity ||23|| . More recently, Farrar proposed |2(J the 
neutral hadron S°, a guds bound-state, as the lightest (/-hadron (see also the calculations 
in the MIT bag model of Ref. ||). 

There is some controversy if a light gluino, with a mass of a few GeV, is allowed. As it 
stands, the Farrar model j^Oj is in conflict with searches for glueballino decays |^5| |26|, [27 



as well as for decays of other unstable g-hadrons PSI. However, these searches had been 



restricted to a narrow band of lifetimes and masses, and their results are not valid in the 
context of more generic models. 

The existence of a light gluino {m g < 5 GeV) can be (dis-) proved due to its contribu- 
tion to the running of a s and to QCD colour coefficients in a practically model-independent 
way. The authors of Ref. |29| used the ratio R between the hadronic and the jl/i produc- 



tion cross-section in e + e annihilation at different energies to constrain the light gluino 
scenario. They excluded light gluinos with mass rrig = 3(5) GeV with 93(91)% CL, while 
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the mass range < 1.5 GeV remained essentially unconstrained. Combining these results 
with the determination of QCD colour coefficients from the analysis of multi-jet events 
30| , the conclusions of pS|| became much stronger: light gluinos with mass < 5 GeV 



in 



were excluded with at least 99.89% CL. The analysis of multi-jet events relied however 
on the use of Monte Carlo simulations which parameters are tuned to QCD without light 
gluinos. Moreover, the multi-jet analysis was based on a tree-level calculation with rather 
large scale ambiguities. The assessment of these uncertainties is difficult, thus preventing 
the definite exclusion of a very light gluino by this argument [El], j3T 



Direct accelerator limits for the gluino as LSP were discussed recently in Refs. [22, p 



The authors of Ref. |33|] concluded that the range 3 GeV < m g < 130-150 GeV can be 
excluded at 95% C.L. based on currently available OPAL and CDF data. Their results are 
sensitive to the details of the hadronic interactions of g-hadrons and, for certain choices 
of the parameters, a window in the intermediate mass region 23 GeV < m g < 50 GeV 
remains open. Meanwhile, Ref. [^] noted that these limits could be weakened if squarks 
are not very heavy and contribute to the jet + missing energy signal, while Ref. f34] 
confirmed an open window for a gluino with 25 GeV < m g < 35 GeV. 

We also mention here that cosmological constraints do not exclude both light and 
heavy gluinos of interest [|17j, |3^] . 

The Gustafson experiment |57j] does not exclude g-hadrons (see section 5). 

Till now we discussed the limits on the gluino mass m g . The lightest g-hadron with 
mass M g is heavier than the gluino by the mass of its constituent gluon or quarks, which 
are expected to be less than 1 GeV. 

Light g-hadrons with M g ~ 1.5 GeV have a spectrum with the GZK cutoff beyond the 
currently observed energy range (see [17] and Fig. of the present paper). Together with 
the accelerator limits on gluino masses this leaves a narrow band of allowed masses for the 
light g-hadrons at 1.5 < M g < 4 GeV. But this window is closed if, as argued in [23|, the 



charged gluebarino guud is lighter than the neutral guds. Indeed, production of charged 
gluebarinos in the Earth atmosphere by cosmic rays and their accumulation in the ocean 
results in too high abundance of "wild hydrogen" in contradiction with observational 



data. In Refs. |2(J, |24|, however, it is argued that the lightest gluebarino is the neutral 
flavour singlet guds, due to strong quark attraction in this state. But even in this case 
the restriction J23| might work, if gwds-gluebarino and proton are bound into anomalous 
deuterium. 

In conclusion, a light gluino — although being disfavoured by various arguments — is 
not excluded. We are studying here the interactions of gluinos, being inspired by a 
possible correlation between AGN and the arrival directions of UHECR j|, and by 
the recent suggestion [EJ, E_8|, that Extensive Air Showers (EAS) observed at the highest 
energies could be produced by g-hadrons. The authors of Ref. [Ej| performed a detailed 
simulation of EAS induced by g-hadrons, using however a phenomenological description 
for the g-hadron-nucleon interaction which is not self-consistent. They found that masses 
as high as 50 GeV are compatible with presently available data. In contrast, it was argued 
in Ref. [17], using kinematical arguments, that the observed shower characteristics exclude 
any strongly interacting particle much heavier than a few GeV. 

The purpose of the present work is to study EAS produced by g-hadrons and to restrict 
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the mass range in which the g-hadron is a viable UHE primary using a self-consistent 
interaction model. We have used for the simulation of air showers initiated by g-hadrons 
a suitably modified version of the QGSJET model |39], [40] which is known to describe 
successfully proton air showers |39], [4l[| . Specifically, we have considered the glueballino 
as the g-hadron but we expect that our results apply to all g-hadrons. We paid special 
attention to consistent calculations of glueballino-hadron interaction cross-sections and of 
cascade particle production in the atmosphere. We have found that the development of 
showers initiated by g-hadrons with masses above 5 GeV differs substantially from proton- 
initiated showers and is inconsistent with the current experimental data. In the window of 
masses 1.5-4 GeV, where a g-hadron can be allowed, glueballino-induced showers do not 
contradict the available data. Future observations by the detectors HiRes and Auger can 
either confirm or exclude g-hadrons as the dominant primary combining the information 
from shower profiles and energy spectrum. However, the best way to test this hypothesis 
is a modified Gustafson experiment (see Section 5). In the case of the discovery of light g- 
hadrons in such an experiment we shall reliably know their properties, thus enabling us to 
calculate the production of these particles in astrophysical sources and their detection in 
the Earth atmosphere. Their absence will preclude further discussion of this hypothesis. 



2 Glueballino— proton (nucleus) interaction 
2.1 QGSJET framework 

QGSJET, a Monte Carlo generator of hadron-hadron, hadron-nucleus, and nucleus-nucleus 
interactions [^, [l(J, was developed in the framework of the Gribov-Regge approach and 
is based on the quark-gluon string model of the supercritical Pomeron |[42|| . Hadronic 
interactions are described as a superposition of elementary rescattering processes between 
the partonic constituents of the projectile and target nucleons (hadrons), resulting in the 
production of colour neutral strings, which further fragment into secondary hadrons. The 
key parameters of the model are the intercepts and the slopes of the Regge trajectories 
of the Pomeron and of secondary Reggeons. These parameters govern the formation of 
different interaction configurations, how the energy-momentum is shared in elementary 
interactions, and also the string hadronization. The model was generalized to hadron- 
nucleus and nucleus-nucleus interactions in the framework of the Glauber-Gribov ap- 
proach E3|, |44J , taking into account low mass diffraction and inelastic screening processes 
|45| . Hard QCD processes were included into the Gribov-Regge formalism via the concept 



of a "semihard Pomeron", which is a t-channel iteration of the soft Pomeron and a QCD 



parton ladder contribution [40, 46 



QGSJET describes hadron-hadron interaction amplitudes as a sum of two contribu- 
tions, namely soft and semihard rescatterings ||40| . The soft contributions are of purely 
nonperturbative nature and correspond to the case of a parton cascade with virtualities 
smaller than some cutoff Qq. Below this cutoff, perturbative QCD is not applicable and 
the interaction is described by phenomenological soft Pomeron exchange. The amplitude 
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fa C for Pomeron exchange between two hadrons a and c is given by [[42 



fl(s,b) = ( 1)2 



KM = R 2 a + R 2 c +a' p y, (2) 

where y = In s is the rapidity size of the Pomeron, s is the squared center-of-mass energy 
for the interaction, b is the impact parameter between the two hadrons, and the parameters 
7a(c)> -^a(c) are the vertices and slopes for the Pomeron-hadron a(c) coupling, respectively. 
Finally, A and aL are the parameters describing the overcriticality and the slope of the 
soft Pomeron trajectory. 

Contrary to soft rescatterings, semi-hard ones correspond to the case when at least 
a part of the parton cascade develops in the region of parton virtualities q 2 > Qq and, 
therefore, can be described on the basis of QCD techniques. The complete semi-hard con- 
tribution is represented by a QCD parton ladder sandwiched between two soft Pomerons 
f46| . For the Pomeron, the formulas can still be used. However, since it is now 

coupled to a hadron a(c) on one side but to a parton ladder on the other side, the slope 
Rc( a ) an d the coupling 7 c(a ) have to be replaced by the slope Rf^ and the coupling V y ^ d of 
the Pomeron- ladder. The latter is parameterised as 

vZ(y) = r[l-ex I (-y)f P , (3) 

where the parameters r and /3p describe the momentum distribution of a parton (sea 
quark or gluon) in the soft Pomeron. Using R 2 ad ~ 1/Qq <C R 2 a (b) + a pV> the slope -Rj 2 ad 
can be neglected. 

To complete the description of soft and semihard contributions, the momentum distri- 
bution function for soft Pomeron emission by a hadron of type a has to be specified. 
It is parameterised in the form 

where the first factor does not depend on the hadron type and describes the probability 
to slow down the hadron constituents to which the Pomeron is connected. In QGSJET, 
the Pomeron is connected to a (dressed) quark-antiquark pair. Using Regge asymptotics 



42| , i.e. a g q ~ 0.5 as the intercept of the Regge gg-trajectory for light quarks, it follows 

a = 1 — 2a qq ~ . (5) 

Similar, the second factor in Eq. (f|) describes the probability to slow down the "leading" 
hadron state configuration; the parameter (3 a is expressed via the intercept a aa of the 
corresponding Regge trajectory as (3 a = —a a a- 

The semihard contribution described above corresponds to the case that gluons or sea 
quarks of the hadron start the interaction at the initial scale Qq. Additionally, valence 
quarks can interact with q 2 > Qq. Then the only nonperturbative input needed are the 
valence quark momentum distributions q v (x, Qq) at the initial scale Q^. 
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2.2 Extension for glueballino 

The nucleon-glueballino interaction can be treated in the QGSJET model in the same 
framework as the one for usual hadrons |39| . The main difficulty is to connect the unknown 
physical parameters (coupling 73, slope and momentum distribution Ng) describing 
the interactions of glueballinos with the corresponding quantities of usual hadrons. We 
use simple scaling arguments to derive the glueballino parameters from those of the pion: 

1. The coupling 7 a of a hadron a to the Pomeron depends essentially on its size and, 
consequently, on its reduced mass M a . If r a denotes the radius of the hadron a 
with the reduced mass M a , then 7 a ~ r\ ~ M~ 2 , where we have neglected a factor 
a s (M 2 ) in the last step. Thus, the Pomeron-glueballino vertex 7^ can be expressed 
via the Pomeron-pion vertex j n as 

76 = 7 -(l|) • (6) 
For the reduced mass of the glueballino, we use 

rrig + m g 

where m g is the mass of the gluino and m g ~ 0.7 GeV is the constituent mass 
of the gluon. Similar, we use for the pion M n = m q /2 with m q ~ 0.35 GeV as 
quark constituent mass. Note that Eq. (|) does not take into account the different 
colour factors of quarks and gluon/gluinos because we consider an effective Pomeron 
coupling to the hadron as a whole, not to individual parton constituents. 

2. The slope R 2 ^ for the Pomeron-glueballino coupling is also inverse proportional to 
M?. Therefore, R 2 ^ is small compared to R 2 and Rp = a'py and can be neglected 
in the formulas (|l]-0), 

XaM^R-l + a'py. (8) 

3. The momentum distribution for Pomeron emission is again given by Eq. (^), Ng ~ 
(xp) a (l — Xp)^o. Now the "leading" configuration consists of a valence gluon and 
gluino, 

P6 = l + Pg + /3 g . (9) 

Assuming that a valence gluon behaves similar to a valence gg-pair in the low x-limit 
gives (3 g ~ 1 — 2a q q ~ 0. The remaining unknown parameter j3 g can be found from 
the momentum distribution between the valence constituents of the glueballino. 
Using as ansatz for the momentum distribution pj? of the gluino 

Pi M ~ 4" (1 - •'',)''" (10) 

and assuming that the energy is shared according to the constituent masses of the 
valence partons, we obtain for the average momentum fraction carried by the gluino 

/ \ m 9 /% + 1 fft\ 

mg + m-g fig + f3g + 2 
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This results in 



/% = ^(& + l)-l^-l. (12) 



m g m g 



Having fixed the free parameters describing the Pomeron-glueballino interactions us- 
ing essentially only one simple, physically well- motivated scaling argument, the soft and 
semi-hard contributions are determined. These two contributions to the total nucleon- 
glueballino interaction are referred below as the contribution due to the soft coupling, 
because they are both caused by soft Pomeron emission of the glueballino. 

To complete the formalism, we need to define the momentum distributions of the 
valence gluon or gluino inside the glueballino probed at the initial scale Qq, when they 
are involved into hard interactions. We shall refer to this contribution below as the the 
contribution due to the direct coupling. Parton emission by a gluino in the s-channel is 
strongly suppressed kinematically in the nonperturbative region q 2 < Qq by its mass: the 



virtuality q 2 of the process g —>■ g + g is determined by the off-shellness q\ — M~ 
produced t-channel gluino, 



of the 



,2 



q 2 -M 2 =^- + M 2 (l-z). (13) 



Here p± and z are the transverse momentum and the light cone momentum fraction for the 
t-channel gluino. Therefore, the gluino momentum distribution at the scale Q\ essentially 
remains in the form ([TOf) , 

9v(x- g ,Ql) =pf(x s ) . (14) 



By contrast, the momentum distribution of the valence gluon differs from p? (1 — x g ) in 



fllPf) because of the emission of sea quarks and gluons (due to the soft coupling defined 
above), 

g v (x g ,Ql)^x^(l-x g )^ +s . (15) 

The parameter 5 is fixed by the momentum conservation constraint for the complete 
(valence and sea) parton momentum distributions at the initial scale Qq. This completes 
the formulation of the initial conditions for the perturbative evolution. 

The treatment of the perturbative part of the interaction is performed to leading-log 
accuracy, the corresponding techniques are described in Ref . |$9[ . All emitted (s-channel) 



partons undergo timelike cascading according to the standard algorithm with soft 
gluon coherence taken into account. At the final stage, soft strings are assumed to be 
formed between on-shell partons according to their colour connection pattern. The final 
gluino is assumed to pick up a gluon-gluon singlet pair from the vacuum. After a colour 
rearrangement similar to the one in Jjip production, a glueballino is formed. String 
hadronization completes the procedure for glueballino-hadron interaction. The extension 
of the model to glueballino-nucleus collisions is based on the standard Glauber-Gribov 



approach and does not differ from the usual hadron-nucleus case [|S . 
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2.3 Numerical results 



The model developed in the last subsection allows both to calculate the cross-sections 
for glueballino-nucleon interactions and to treat consistently particle production in these 
reactions. Some quantities characterising the glueballino-nucleon interactions are given 
in Table 1 for E Uh = 100 GeV and in Table 2 for £ lab = 10 12 GeV. We present both 
total and inelastic cross-sections as well as the partial contributions arising due to the 
soft (cr^ coupL ) and direct (o"f ^ coup1 ' ) coupling of the glueballincQ. At low energies, the 
soft coupling strongly dominates the G-proton interaction for all M G being governed by 
nonperturbative soft interactions, while the direct coupling can be neglected. At high 
energies, this picture changes considerably. The soft coupling becomes more and more 
suppressed for large M G . By contrast, the direct contribution, which is purely perturbative 
on the glueballino side, is nearly independent of M G . This important difference from 
the usual hadron case is due to the very asymmetric energy partition between parton 
constituents of the glueballino. For large M G , the valence gluino carries almost the whole 
initial energy of the particle (88% for M~ g = 5 GeV and 99% for M- g = 50 GeV) - Eqs. 
(|T0] - |TT|) , flH| ), thus leaving just a small part of it to other partons, to which Pomerons are 



connected. Therefore, the glueballino behaves in the limit of large M G essentially as a 
perturbative object, as one expects from kinematical considerations |T7[. Finally, the last 



raw of the tables shows the inelasticity coefficient K me \ as another important quantity 
which distinguishes proton-proton and G-proton interactions. Although at energies of 
interest for UHECR the total cross-section for G-proton interactions is rather large, a 
heavy glueballino behaves like a penetrating particle in the atmosphere loosing only a 
small part of its energy in one interaction. This conclusion was already reached in Ref. |17] 



from semiqualitative considerations. The reason for this effect is twofold. On one side, 
as discussed above, gluinos of larger masses carry a larger fraction of the initial particle 
energy, leaving a smaller part of it for the sea constituents ((anti-) quarks and gluons) and 
thus reducing the average number of multiple interactions in G-proton (nucleus) collisions. 
On the other hand, the relative weight of the "direct hard" process increases for heavier 
gluinos, where the valence gluino loses typically only a small part of its energy, as a 
large longitudinal momentum transfer is strongly suppressed in that case by the process 
virtuality q 2 



^ z 2 M? 



«^ = i^ + tf*. (16) 

with p± and z being the transverse momentum and the light cone momentum fraction for 
the t-channel gluon emitted of the initial valence gluino which mediates the gluino hard 
interaction with the target proton (nucleus). 

We show also in Figs. [I] and ^| the total and inelastic cross-sections as function of the 
interaction energy for glueballino masses M G = 2,5 and 50 GeV. The fast increase of 
the direct contribution with i? lab produces an interesting effect: the total interaction cross- 
section for the largest glueballino mass considered, M G = 50 GeV, which is dominated by 



1 Note that the total cross section is smaller than the sum of the partial contributions (7t ot coupL and 
°tot COUpl ! because it contains contributions from multiple rescatterings of both types. 
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the direct contribution, overshoots the ones for smaller glueballino masses in the energy 
range 10 4 - 10 9 GeV. 



3 Extensive air showers (EAS) 

In this section, we present some results of our simulations for the glueballino-initiated 
EAS. Shower profiles and distributions of shower maxima are shown in Figs. |3|-|8]for three 
different initial energies, Eq = 10 17 , 10 19 and 10 20 eV, for the glueballino as primary with 
different choices of the gluino mass. For comparison, the case of a primary proton is also 
shown. At the highest energy considered, E = 10 20 eV, the longitudinal shower profiles 
(Fig. H) and the distribution of the shower maxima X max (Fig. |) of glueballino-induced 
EAS are comparable with those induced by protons in the case of glueballino masses 
smaller than 5 GeV. As the glueballino mass increases, the shower develops deeper in the 
atmosphere with a less pronounced maximum. The fluctuations in X" max increase also for 
larger glueballino masses. The main reason for both effects is the competition between 
the large glueballino-nucleus cross-section and the small inelasticity of the interactions: a 
heavy glueballino injects in one interaction only a small part of its energy into secondary 
hadronic and electromagnetic cascades, while interactions with a large momentum transfer 
are rare and increase only the fluctuations. Figure [| clearly shows that the glueballino- 
induced EAS drastically differ from the proton-induced showers for glueballino masses 
larger than 5 GeV, and hence these showers can be distinguished even in case of low 
statistics. In case of a lighter glueballino with M G = 2 GeV, a larger statistics is necessary 
to distinguish glueballino from proton, when only X max measurements are used. The same 
conclusions can be drawn comparing the shape of the calculated profiles for individual p- 
and G-induced EAS of energy E = 3.2 • 10 20 eV with the corresponding measurements of 
the Fly's Eye collaboration j|8], cf. Fig. |9]. In doing so we choose only those showers which 



reach their maxima near the measured value X max = 815±50 g/cm 2 . Then we average the 
obtained profiles and shift them to the same position of the shower maximum, X max = 815 
g/ cm 2 . It is easy to see that for gluino masses larger than 5 GeV the shape of the calculated 
profile strongly disagrees with the experimental observations. The account for the LPM 
effect results in only 5% reduction of the electron number in the shower maximum for 
proton-induced EAS [19j and has an even smaller influence on the G-induced showers due 
to much softer 7r°-spectrum in the glueballino interactions. 

The calculated lateral distribution functions (LDFs) for electrons and muons (E^ > 1 
GeV) at the AKENO observation level (900 g/cm 2 ) are shown in Figs. [TTJ and |ll| for the 



proton and glueballinos with masses 2 and 5 GeV. The plotted values are the LDFs of 
electrons and muons p e {R), P^{R) &t different distances R from the shower core. Although 
these distributions are substantially different for showers initiated by glueballinos and 
protons, they hardly can be used to search for the light glueballino on the basis of existing 
data, e.g., of AGASA. An adequate tool for the glueballino search is the fluctuation of the 
muon density at distances R > 300 m from the core. This method allows in principle to 
discriminate showers initiated even by light glueballinos with Mq « 2 GeV from proton- 
initiated EAS (see Table |j). 
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Finally, we shall compare our results with those of Albuquerque, Farrar and Kolb 
(AFK) p8| . AFK have modified the event generator SIBYLL including the g-hadron 
(G) as a new particle. The interaction properties of g-hadron were taken ad hoc. Two 
assumptions were used for the total cross-section: <J tot (Gp) ~ cr tot (irp) (the favourite 
choice) and cr to t(G'p) ~ 0.1o"tot(7rp)- The mean energy fraction transferred from the g- 
hadron to the shower per interaction was modeled by a Peterson fragmentation function. 
Hard interactions with the production of minijets by the incident g-hadron were neglected. 

It is easy to see that these modifications are not self-consistent. Indeed, on one hand 
the authors assume a large Gp cross-section, while on the other hand they neglect the 
hard processes (production of minijets), which give the dominant contribution to the 
Gp cross-section and make it large. In fact, our calculations explicitly show that at 
ultra-high energies the Gp cross-section and particle production are dominated by hard 
interactions for both light and heavy gluinos. For light gluinos, valence gluon and sea 
partons have enough momenta for hard interactions. For heavy gluinos, its own ("direct") 
hard interaction dominates. Soft interactions without the production of parton jets are 
negligible in both cases. 

To elucidate the reason for the failure of the AFK approach, we have calculated the 
total G-nucleon cross-section switching off the hard interaction (cf. 5th entry, oafk; in 
Table 1 and 2): At energies relevant for UHECR, the interactions considered by AFK are 
only subdominant and result in much smaller total cross-sections as compared with ours 
or those assumed by AFK. 

4 Energy losses of glueballinos on CMBR photons 
and glueballino energy spectrum 

Although both valence constituents of the glueballino are electrically neutral, UHE glue- 
ballinos loose energy due to scattering on CMBR photons. The value of the cutoff in its 
energy spectrum is determined by the transition from adiabatic energy losses (redshift) 
to rapidly increasing energy losses due to the reaction G + 7 — > G + 7r° at higher en- 
ergies. This process cannot occur due to 7r° exchange in the t-channel. The dominant 
contribution is given by the resonant formation of gqq states in the s-channel. The mass 



spectrum of gqq states was calculated as function of the gluino mass in Ref. |49| in the 
MIT bag model. The lowest gqq state found was the spin- 1/2 state p\j%\ its mass difference 
to the glueballino is however, except for itlq < 1.2 GeV, too small as to allow the decay 
P1/2 — ► G + tt°, cf. Table |6|. We assume therefore that the first resonance in the s-channel is 
an excited py 2 state; for its mass m(py 2 ) we use m(py 2 ) = m(pi/ 2 ) + 730 MeV guided by 
the mass difference between the p(770) and the p(1400). The Breit-Wigner cross-section 
for the reaction G + 7 — > px/2 — ► G + ir° is p2| 



/ x _ 2tt B in B out Tf ot mj 

where p cm and yfs are the momentum and the total energy of the particles in the cm 
system, m~ p is the mass of the P1/2/ 1 P1/2 particle and r tot is its total decay width. Finally, 
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B m and B out denote the branching ratios of the resonance to the initial and final states. 

For the determination of the unknown branching ratios B m , B ont and the total decay 
rate r tot , we can use the analogy between the glueballino and the pion together with 
scaling arguments as in Sec. |^ both for the py 2 and the py 2 . Then 



B out = T{~ Pl/2 + 7r°)/r tot (p 1/2 ) ~ 1 



and 



-Bit 



1/2 



G + i) r( P ° 



7T 



+ 7) 



Ttot (P1/2) 



r to t(p) 



7- 10" 



(19) 



The total decay rate r tot (pi/ 2 ) 
with 



oc gimp can be calculated using rrip from Ref. |49] together 

(20) 



9p 



m; 



9 P - 



The maximum of the photo-pion production cross-section for the glueballino is at least 
a factor 8 smaller than for the proton due to the differences in B[ n , the spin factors and 
Pcm- We neglect multipion production G + 7 — > G + n + + n~ , which operates at energies 
above those we are interested in. The non-resonant contributions have been calculated in 
the vector dominance model and give only a small contribution to the total cross-section 
of the glueballino. The hard processes with gluon exchange are important only at high 
energies. 

The energy loss of a particle scattering on CMBR photons is given by |]50 



1 dE 
E ~dt 



T 



27r 2 r 2 J Eth 



dE r a(E r )y(E r )E T 



1 — cxp 



E r \ 
2TTJ 



(21) 



Here, T is the Lorentz factor of the particle in the CMBR frame, E r is the energy of 
a CMBR photon in the rest system of the glueballino, _E t h = m^il + m^jlm^) is the 
threshold energy in the glueballino rest system and y is the average fraction of energy lost 
by the glueballino in a collision, which for one-pion production is 



y(E r 



E r 1 + ml/2E r rriQ 
tjiq 1 + 2E r /mQ 



(22) 



The energy spectrum dN/ dE of glueballinos emitted by diffuse sources is at the present 
redshift z = given by 



dz a (1 + z a 



-5/2 



dN dE g (E,z g 
dE f } 9) dE 



(23) 



where dN/dE g (t g ) is the injection spectrum at redshift z g and a is a constant depending 
mainly on the source emissivity. We have calculated dE g (E, z g )/dE as it is described in 
Ref. ||51]| . For the injection spectrum we have used dN/dE g (z g ) ~ E g 2J with a maximal 
redshift z max = 2 and both without an intrinsic energy cutoff (Fig. |12]) and with a cutoff 
at E = 1 • 10 22 eV (Fig. 0). 
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The calculated diffuse spectra are presented for three different masses of the glue- 
balhno, m G = 1.5, 2 and 3 GeV, and for comparison the corresponding proton spectra 
are also shown. The cutoff in the gluebalhno spectra is shifted to larger energies not 
only because of m G > m N but also because of the large mass gap between G and p\, 2 . 
Moreover, the smaller cross-section and energy transfer make the cutoff less pronounced. 
The calculated spectra are for all three masses in agreement with the observations |36| . 
However, the observed flux above the GZK-cutoff can only be reproduced, if either the 
gluebalhno injection spectrum is rather flat or if also at energies E ~ 10 19 eV glueballinos 
are a non-negligible component of the cosmic ray flux. Finally, we stress that the exact 
form of the gluebalhno energy spectrum depend rather strongly on the mass spectrum of 
the hadronic bound states of the gluino; still our general observations (shift of the cut- 
off to larger energies, less pronounced cutoff) are independent of the details of the mass 
spectrum. 



5 Experimental limits from beam-dump experiments 

The experiment FNAL-E-0330 by Gustafson et al. was designed to search for neutral 
hadrons with masses > 2 GeV and lifetimes > 10~ 7 s 113711 . In contrast to the recent 



experiments |25|, |26], [27], [2q| , the Gustafson experiment was sensitive to long-lived or 
stable hadrons. It set upper limits on the production cross-section of these hadrons in the 
mass range 2-12 GeV which are given in the second line of Table [|. However, the relation 
between the measured mass, M meas , in the Gustafson experiment and the physical mass 
of g-hadrons needs careful consideration. 

In the Gustafson experiment, short pulses of protons with energy 300 GeV produced in 
the beam-dump target neutral secondary particles, which were detected in a calorimeter 
at distance I = 590 m from the target. The calorimeter had the total thickness X ca \ = 
900 g/cm 2 , and it was assumed that a particle looses all its energy therein. Then the 
measured difference in the time of flight 5t of massive secondaries and photons determines 
the mass of the secondary, 

M meas = (2cElJt/l) 1 / 2 . (24) 

Here, the energy loss Ei oss in the calorimeter is taken as the energy of a particle. In 
fact, a g-hadron looses however only a part of its energy in the calorimeter: E\ oss /E = 
K inel X cld /X int , where X int is the inelastic interaction length. Then the physical mass M G 
of a g-hadron is connected with the measured mass as 

M meas = (E loss /E)M G . (25) 

In the case of a gluebalhno, one obtains from Eq. (^) together with the data of Table [l| 
^meas ~ 1 GeV for M G = 2 and 5 GeV. Thus, in both cases glueballinos fall in the region 
of neutron background and therefore these masses are allowed. 

Are heavier masses of g-hadrons allowed by the Gustafson experiment? 

We have calculated the cross-section for gluino production in pp-collision via the two 
subprocesses gg — > gg and qq — > gg at next-to-leading order |52| for £i ab = 300 GeV 



using the program Prospero [53|. As parton distributions we have chosen those of Gliick, 
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Reya and Vogt [|54[], while we have used Q = m g as renormalization and factorization 
scales. In Table |5|, the production cross-sections averaged over the region allowed by the 
experimental cuts are compared with the upper limits from the Gustafson experiment. 
The inspection of Table [| for larger masses shows that the calculated cross-sections are 
smaller than the upper limits given in the second line. We conclude that the results of 
the Gustafson experiment do not exclude g-hadrons. 

A modified Gustafson experiment has great potential to discover light g-hadrons in the 
window 2-4 GeV or to reliably exclude it. For this aim, a thicker calorimeter is needed. 
If g-hadrons lose all their energy in the calorimeter, then M meas « M G and the third 
and fourth line of Table [| show that the calculated cross-sections for Mq in the interval 
2-4 GeV are larger than the experimental upper bounds. It means that g-hadron in this 
mass window can be discovered or excluded, if the calorimeter is thick enough. Taking 
into account the large penetrating power of g-hadrons the neutron background can be 
greatly reduced by placing the absorber behind the target which thickness is tuned to 
absorb the neutrons but to be transparent for g-hadrons. 

The signature of g-hadrons is given by the compatibility of the gluino production 
cross-section with the measured flux of detected particles. The path-length and the av- 
erage fraction of lost energy will serve as further indicators. The measured properties 
of g-hadrons will then allow reliable calculations of g-hadron production in astrophysical 
sources and their detection in EAS. 



6 Conclusions 

There are two mass windows where a light gluino is marginally allowed: m g < 3 GeV 
and 25 < m g < 35 GeV. The first window is disfavoured by the gluino contribution to 
the running of a s and to the colour coefficients. The lightest g-hadron is heavier than 
the gluino by the mass of constituent gluon or quarks. The existence of a light quasi- 
stable g-hadron, corresponding to the first window for gluino masses, crucially depends on 
whether the neutral or charged g-hadron is the lightest one. A light quasi-stable charged 
g-hadron is forbidden by the production of "wild hydrogen" in the Earth atmosphere by 
cosmic rays. The case that the lightest g-hadron is neutral is also forbidden, if it forms 
a bound state with the proton ("wild deuteron"). The status of the light gluino and 
lightest g-hadron can further be clarified by further theoretical analysis. At present we 
conservatively consider the first gluino window as disfavoured, but not excluded. 

As the carrier of UHE signal, the light g-hadrons with a mass larger than 1.5 GeV have 



a spectrum with the GZK cutoff beyond the observed energy range (see Fig. 12). Together 
with accelerator limits on the gluino mass, m g < 3 GeV, it leaves for g-hadron masses 
the narrow window 1.5-4 GeV. The other window allowed by accelerator experiments is 
25-35 GeV. 

In this paper we have studied the interaction of g-hadrons with nucleons and nuclei 
and the development of UHE EAS initiated by such particles. In practice, we have 
considered the special case of a glueballino as the primary particle. We think, however, 
that any other g-hadron with equal mass has essentially the same interaction properties. 
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We have calculated the glueballino-nucleon inelastic and total cross-sections for different 
masses of glueballinos (see Figs, [l], H). The longitudinal development of the EAS with 
energy 3 • 10 20 eV in the atmosphere is shown in Fig. [| One can see that glueballinos 
heavier than 5 GeV resemble penetrating particles. The profiles shown in the Figure 
are directly measured in Fly's Eye experiment, and the data are quite different from the 
graphs displayed here for heavy glueballinos. The second mass window 25-35 GeV can be 
already reliably excluded on the basis of these measurements, though a detailed analysis 
is desirable. 

The showers produced by g-hadron from the low- mass window are similar to proton- 
initiated showers in all characteristics (see the longitudinal profiles in Figs. |3], |5| and ^, the 
distributions over X max in Figs. f|, || and ||, lateral distributions for electrons and muons 
in Figs. [10], |TT] and fluctuations in muon and electron densities in Tables || and In 



principle, the best possibility to distinguish the showers produced by g-hadron with mass 
2 GeV from the proton-induced showers are the fluctuations in the muon density at large 
distances d > 600 m (see Table |3|). However, the statistics of the largest array at present, 
AGASA, is not sufficient for such a discrimination. The large A max tail in the distribution 
of showers initiated by g-hadrons over X max offers for HiRes or AUGER another, more 
promising possibility to (dis-) prove light gluinos as UHE primaries. 

More reliably the quasi-stable 2-3 GeV g-hadron can be found in an especially designed 
accelerator beam-dump experiment (see section 5). In this experiment, it is possible to 
discover supersymmetry (light gluino) and to find a new primary for the UHE cosmic 
signal. 
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7T 


M G = 2 GeV 


M G = 5 GeV 


Mq = 10 GeV 


M G = 50 GeV 


Ctot 

0"in 

s— coupl 
^tot 

d— coupl 
^tot 

0"AFK 
^inel 


24.4 
20.9 


3.9 
3.8 
3.8 
0.14 
3.8 
0.25 


2.9 
2.8 
2.7 
0.15 
2.7 
0.15 


2.4 
2.4 
2.4 

0.013 
2.4 

0.074 


2.2 
2.1 
2.2 


2.2 
0.0054 



Table 1: Total cross-section cr tot , inelastic cross-section a m , cross-section due to soft 
coupling <7to^ coupI , cross-section due to direct coupling cr t d ^ coupl , cross-section without hard 
interactions o"afk and the inelasticity coeffiecent K me \ for the glueballino-nucleon scatter- 
ing. The data are given for four values of the glueballino mass M G from 2 GeV to 50 GeV 
at glueballino energy E\ a ^ = 100 GeV. As references, the total and inelastic cross-section 
of the pion are also given. All cross-sections are in mbarn. 





7T 


M G = 2 GeV 


M G = 5 GeV 


M G = 10 GeV 


M G = 50 GeV 


Ctot 


152 


103 


94.9 


91.3 


88.0 


0"in 


112 


71.3 


65.9 


63.4 


60.2 


s— coupl 
Otot 




85.5 


69.4 


61.8 


50.8 


d— coupl 
^tot 




72.1 


72.4 


72.3 


74.4 


0"AKF 


92.2 


18.2 


13.3 


11.8 


10.6 


^inel 




0.26 


0.14 


0.082 


0.018 



Table 2: The same as Table 1 for E ]ah = 10 12 GeV. 



R, m 


10 


50 


100 


200 


300 m 


600 m 


1200 m 


proton 


0.077 


0.054 


0.068 


0.095 


0.11 


0.15 


0.18 


G(2 GeV) 


0.20 


0.25 


0.27 


0.30 


0.31 


0.34 


0.36 


(7(5 GeV) 


0.33 


0.37 


0.39 


0.41 


0.43 


0.45 


0.47 


G(10 GeV) 


0.40 


0.43 


0.45 


0.47 


0.48 


0.50 


0.52 



Table 3: Normalized la fluctuations of the electron LDF a Pe (R)/ p e (R) at the AKENO 
observation level for EAS of energy Eq = 10 20 eV initiated by protons and glueballinos of 
masses 2, 5 and 10 GeV. 
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R, m 


10 


100 


200 


300 m 


600 m 


1200 m 


proton 


0.19 


0.16 


0.16 


0.17 


0.20 


0.25 


G(2 GeV) 


0.18 


0.20 


0.24 


0.27 


0.33 


0.41 


G(5 GeV) 


0.23 


0.30 


0.34 


0.37 


0.43 


0.51 


G(10 GeV) 


0.30 


0.36 


0.40 


0.42 


0.48 


0.55 



Table 4: Normalized la fluctuations of the muon LDF (E^ > 1 GeV) v Ptl (R)/ Pn{R) at 
the AKENO observation level for EAS of energy Eq = 10 20 eV initiated by protons and 
glueballinos of masses 2, 5 and 10 GeV. 



M meas /GeV 


2 


3 


4 


6 


8 


10 


90% C.L. 


1.0e-32 


3.6e-33 


1.3e-33 


5.3e-33 


2.0e-33 


1.7e-33 


rrig/GeV 


2 


3 


4 


6 


8 


10 


Ed 3 a/d 3 p 


3.8e-30 


1.7e-31 


l.le-32 


4.2e-35 


3.9e-38 


8.9e-43 



Table 5: The upper 90% C.L. limit on the invariant cross-sections Ed 3 a/d 3 p (second line) 
for measured masses M meas (first line) given by the Gustafson experiment. The fourth 
line gives the calculated invariant cross-sections for gluino production in pp-interaction at 
E p ,\ah = 300 GeV and for gluino mass m g (third line). All cross-sections are in cm 2 /GeV 2 . 



nig /GeV 


0.64 


1.04 


1.48 


2.41 


M d /GeV 


1.00 


1.50 


2.00 


3.00 


m Pi /2 / GeV 


1.17 


1.58 


2.04 


3.03 



Table 6: The masses of the glueballino G and the pi/2 as function of the gluino mass m g 



according to Ref. [§|| . 
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Figure 1: Total glueballino- nucleoli cross-section <x tot in mbarn as function of E\ ah for 
M d = 2, 5 and 50 GeV. 
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Figure 2: Inelastic glueballino- nucleon cross-section <r in in mbarn as function of E\ ah for 
M d = 2, 5 and 50 GeV. 
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Figure 3: Longitudinal shower profile for EAS of energy E 
protons and by glueballinos with M g = 2, 5, 10 and 50 GeV. 
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Figure 4: Normalized distribution p(X max ) of the shower maxima for EAS of energy 
E = 10 20 eV initiated by protons and by glueballinos with M- g = 2 and 5 GeV. 
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Figure 5: Longitudinal shower profile for EAS of energy E 
protons and by glueballinos with M g = 2, 5, 10 and 50 GeV. 
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Figure 6: Normalized distribution p(X max ) of the shower maxima for EAS of energy 
E = 10 19 eV initiated by protons and by glueballinos with M- g = 2 and 5 GeV. 
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Figure 7: Longitudinal shower profile for EAS of energy E 
protons and by glueballinos with M g = 2, 5, 10 and 50 GeV. 
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Figure 8: Normalized distribution p(X max ) of the shower maxima for EAS of energy 
E = 10 17 eV initiated by protons and by glueballinos with M- g = 2 and 5 GeV. 
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Figure 9: Comparison of Flye's Eye highest energy event with the longitudinal shower 
profile for EAS of energy E = 3 x 10 20 eV initiated (from top to down) by protons and 
by glueballinos with M- g = 2, 5, 10 and 50 GeV. The shower profiles are shifted so that 
their X max agrees with the observed shower maximum. 
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Figure 10: Electron LDF p e {R) in m~ 2 at the AKENO observation level as function of 
the distance R (in m) from the shower core for EAS of energy E = 10 20 eV initiated by 
protons and by glueballinos with M g = 2 and 5 GeV. 
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Figure 11: Muon LDF (E^ > 1 GeV) p^R) in m~ 2 at the AKENO observation level as 
function of the distance R (in m) from the shower core for EAS of energy E = 10 20 eV 
initiated by protons and by glueballinos with M g = 2 and 5 GeV. 
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Figure 12: Diffuse glueballino flux from uniformly distributed sources with injection 
spectra dN/dE oc E~ 2 - 7 for Mq = 1.5, 2, and 3 GeV, compared with proton flux and 
observational data from AGASA. Either the gluino flux or a combination of proton and 
gluino fluxes can fit the data. 
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Figure 13: Diffuse glueballino flux from uniformly distributed sources with injection 
spectra dN/dE oc E~ 2 - 7 and intrinsic cutoff E = 10 21 eV for M d = 1.5, 2, and 3 GeV, 
compared with proton flux and observational data from AGASA. Either the gluino flux 
or a combination of proton and gluino fluxes can fit the data. 
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